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ABSTRACT 

We propose a disc-corona model in which a geometrically thin, optically thick disc surrounds 
a Kerr black hole, and magnetic fields exert a time-steady torque on the inner edge of the 
accretion disc. The analytical expression of the total gravitational power is derived from the 
thin-disc dynamics equations by using this new boundary condition. It is shown that the mag- 
netic torque can considerably enhance the amount of energy released in the disc-corona sys- 
tem. Furthermore, the global solutions of this disc-corona system are obtained numerically. 
We find that the fraction of the power dissipated into the corona in the total for such disc- 
corona system increases with the increasing dimensionless black hole spin parameter a*, but 
is insensitive on the Ae which is the additional radiative efficiency parameter relevant to mag- 
netic torque, for As > 1. In addition, the emerged spectra from this disc-corona system are 
simulated by using Monte-Carlo method, and the effect of the different parameters on the 
output spectra is discussed. 
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1 INTRODUCTION 

It is believed that the black hole (BH) accretion disc is an ef- 
fective model for explaining the high energy radiation in astro- 
physics. The standard accre tion disc (SSD) model was proposed by 
Shakura & Sunvaevl dl973l) . in which the disc is geometrically thin 
and optically thick. The SSD is widely used in modeling the spec- 
tral energy distribution (SED) of AGN. The general relativistic SSD 
model has been investigated in detail by Novikov & Throne (1973, 
hereafter NT73), and Page & Throne (1974, hereafter PT74). In 
this SSD model, it has been assumed that there is no stress at disc's 
inner edge, i.e. the so-called "no-torque inner boundary condition". 

Howe ver, the "no- torque inner boundary condition" was ques- 
tioned by Krolikj I I 19991) on the basis th at magnetic fields are th e 
likely agent of the to rque in the discs dBalbus & Hawlevlll998h . 
lAgol & Kroiikl (12000) proposed that magnetic fields connecting the 
disc to the plunging region can exert stresses on the inner edge of 
an accretion disc around a black hole, and they recomputed the 
relativistic corrections to the thin-disc dynamics equations when 
these stresses ta ke the form of a time-stea dy torque on the inner 
edge of the disc. lGammie & Charles also noted that, within 

the confines of a highly-idealized model of inflow dynamics, this 
torque can considerably enhance the amount of energy released in 
the disc. 

Apart from SSD, advection dominated accretio n flow (ADAF) 
is another important model of the accretion flow dNaravan & Yil 
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Il994lll995h . It is believed that black hole X-ray binaries present 
various spectral states, most notably the low/hard state, high/soft 
state and intermediate state. The SSD has been very successful in 
describing the high /soft states of Galactic black hole candidates 
(GBHCs), but their low/hard state characterized by power-law-type 
spectra was problematic to the SSD. The ADAF model is quite suc- 
cessful in reproducing the hard spectra of GBHCs, as well as those 
of low-luminosity AGNs (LLAGNs), however it also has a several 
problems. Thus we cannot expect to interpret the various compo- 
nents of spectra based on one accretion mode. 

The power-law spectra of GBHCs is generally explained by 
Comptonization of softer photo ns by hot electrons fr om a mag- 
netic corona on an accretion disc. lLiang & Price] d 19771) proposed a 
high-temperature corona model analogous to solar corona, and the 
effects on the disc whe n part of the dissipa t ion oc curs in the corona 
were first discussed by 5 onson & Kuperusl ( fl984h . In the accretion 
disc-corona scenario, an optically thin hot thermal corona is located 
above the surface of disc. A fraction of soft photons, which are 
from the cold disc, are Compton up-scattered to X-ray photon by 
hot electrons. In this model the corona can explain the power-law 
X-ray spectra very well, and reprocessing of the coronal X-rays by 
the cold disc gives rise to the observed emission lines naturally. The 
iron K a fluorescence line provides us a diagnostic of the geometry 
of the accretion flow and the property of the space-time around the 
BH. Recent work on the disc-corona model can be fo und in, e.g. 
iMerloni & Fabiart d2002h ; lLiu et alj J2002l) : lcao1 J2008h . 

Shakura & Sunvaevl dl973l) used the famous "a-prescription" 
to deal with viscosity, but the physical process leading to viscos- 
ity and turbulence in the disc remain unclear. The magnetic rota- 
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tion instability (MRI) of weak magnetic fields in accretion discs is 
thought to play an important role in the evolution and dynamics of 
astrophysical accretion discs. lBalbus&Hawievllll998h argued that 
this instability should have a rapid growth rate of the order of the 
orbital frequency Q, resulting in a greatly-enhanced effective vis- 
cosity that is able to transport angular momentum outward. On the 
other hand, the magnetic fields generated in the SSD are strongly 
buoyant, and a fraction of the magnetic energy is transported verti- 
cally to heat the corona above the disc. 

The multi-wavelength observations of nearby LLAGNs have 
revealed that the SEDs of LLAGNs and of GBHCs in the low/hard 
state possess many similarities, e.g. flat, compact radio cores with 
high brightness tempe ratures, and a hard X-ray power-law with 
high a energy cut-off. iMerloni & Fabianl |2002) pointed out that 
there should be a common accretion mode for these low luminosity 
black holes, and proposed a new model for low-luminosity black 
holes, in which a SSD at low accretion rates dissipates a large frac- 
tion of its gravitational energy in a magnetic corona. 

Motivated by the these results, in this paper we shall investi- 
gate a disc-corona model, in which a geometrically thin, optically 
thick disc surrounds a rotating Kerr BH, a magnetic field connect- 
ing the plunging region and the disc exerts a non-zero torque at the 
inner boundary of the disc, and the corona is assumed to be heated 
by the reconnection of the magnetic fields generated by buoyancy 
instability in the disc. This paper is organized as follows. We de- 
scribe the model in section 2, where in a magnetic torque is ex- 
erted on the inner edge of the accretion disc. In section 3 the global 
solutions of this disc-corona system are obtained, and the effects 
of different parameters on the fraction of accretion power dissi- 
pated into the corona are discussed. In section 4 we simulate the 
emerged spectra by using a Monte Carlo method. Finally, in sec- 
tion 5, is the brief discussion. Throughout this paper the geometric 
units G — C = K — 1 are used. 



2 DESCRIPTION OF MODEL 

In our accretion disc-corona system, a geometrically thin and op- 
tically thick disc is sandwiched by a magnetic corona, and part of 
the gravitational energy of the accreted matter is released in the hot 
corona. The general relativistic model for a steady, axisymmetric, 
and thin Keplerian disc around a Kerr black hole has been described 
in detail by NT73 and PT74 . In their model, it has been assumed 
that there is no stress at disc's inner edge. The equation of angular 
momentum conservation is given as follows 
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where M is the accretion rate of the disc,u, ip, p, are the metric 
coefficients, W is integrated shear stress, L + is the specific angular 
momentum of a particle in the disc. 

The total gravitational power dissipated in unit surface area of 
the disc-corona system Q is given by 

Q = {M/4n)e- (v+ ^ +>l) f, (2) 

where the function of radius / is defined by PT74. In the case of 
"no-torque inner boundary condition", the boundary condition on 
/at the radius r ms of marginally stable orbit is f ms = 0. 

If it is taken into account that magnetic stresses exert a torque 
on the inner edge of the accreti on disc, the appropria te boundary 
condition at r ma is expressed as (Agol & Krolik 2000) 
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where Ae is the additional radiative efficiency, Xms = y r ma /M, 
and C is general relativistic correction factor defined by PT74. 

Using this boundary condition we have the new function / of 
radius , that is 



3(E ms -Sl ms L+ ls ) 2 Ae 
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E + and L + are the specific energy and the specific angular mo- 
mentum of a particle in the disc, respectively, and they read 
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X = y r/M is the dimensionless radial coordinate, and a, — 
a/M is the dimensionless black hole spin parameter. 

The power dissipated in the corona (i.e. the magnetic Poynting 
flux in the vertical direction from the thin disc) is 



Q cor — Pm^P 
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where P m is the magnetic pressure in the disc, and up is the ve- 
locity of magnetic flux transported vertically in the disc. Here we 
assume the velocity up of magnetic flux tubes is proportional to 
their internal Alfven velocity, i.e. up = bvA, and b is related to the 
efficiency of buoyant transport of magnetic structure in the vertical 
direction inside t he disc, which is of the or der of unity for extremely 
evacuated tubes jMerloni & Fabianll2002h . 

Now we give the equations of the disc structure as follows. 
The equ ation of vertical pressure b alance in the vertically-averaged 
form is dNovikov & Thornelll973l) 

H = (P/p) 1 ' 2 (r a /M) l/2 AB- 1 C 1/a D- 1 '*E-V a (8) 

where H is the height of the accretion disc, P and p are pressure 
and density of the disc respectively. A,B, C,D, E are general rela- 
tivistic correction factors defined as follows: 

A = 1 + al X ~ A + 1al X -* 
B = l + a tX ~ 3 



D = 1 - 2 X " 2 + al X - A 

E = l + 4a 2 X " 4 - 4a 2 X - 6 + 3atx" 8 



(9) 



The equation of energy conservation is (see Eq.(5.6.13) in 
NT73) 



W = -(M/r 3 )- 1/2 CD' 1 Q, 

and W is integrated shear stress defined as 



W = 2 t rip dz ~ 2t rv H, 
Jo 

where t r u> is the interior viscous stress in the disc. 
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The equation of state for gas on the disc is 

P = P mag + Ptot = P mag + \aT A + p (T/mp) (12) 



where P tot is the total pressure (gas pressure plus radiation pres- 
sure) at disc mid-plane, a is the radiative constant. m v is the rest 
mass of proton, po and T are the density of rest mass and the tem- 
perature in the disc, respectively. 

In the SSD model the interior viscous stress t rlp is usually 
assumed to be proportional to the pressure, i.e. 
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Since the process of generating magnetic fields in the accre- 
tion disc is still unclear, we adopted different magnetic pressures 
as: 



aoPtot - 

OloPgas 
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and as ao is a constant of the order of unity, we adopt ao — 1 in 
our calculations. 

In the disc-corona scenario, the soft photons from the disc are 
scattered by the electrons in the corona to X-ray bands. About half 
of the scattered photons are intercepted by the disc. It is a rough 
approximation because light bending effects ought to be taken into 
account at the vicinity of BH. The reflection albedo ao is relatively 
low (ao ~ 0.1 — 0.2), and most of the inciden t photons from 
the co rona are re-radiated as black body radiation (Zdziars ki et al.l 
1999). Thus the energy transport equation for the disc is 



aT = 2nHp(Q — Q cor + 



ao , 
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where n is the Rosseland mean of total opacity which can be ex- 
pressed as 



k = (0.64 x W^)pT~ 7/2 + 0.40 (cm 2 /g) 



(16) 



and ao = 1.5 is adopted in the calculations. 

Solving equations (2),(7)-(16) numerically, the power dissi- 
pated in the corona Q CO r and the structure of the disc can be derived 
as function of r. The ratio of the power dissipated in the corona to 
the total for such disc-corona system is given as 
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the co ncerned Kerr metric coefficients are given by dThorne et all 
fl986l) 



vj — (S/p) sin 
P 
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3 NUMERICAL RESULTS 

By using Eq.(2)-(6) we have the total gravitational power Q versus 
the radius of the disc for the different values of Ae as shown in 
Fig.l. Note that the power Q is in units of Qo = M/M 2 , where 
M is the BH mass. Inspecting the above Fig, we have the follow- 
ing results: (I) the power Q depends sensitively on the values of 
Ae, i.e. the magnetic torque can considerably enhance the amount 
of energy Q released in the disc; (II) the power Q decreases mono- 
tonically with the increasing radius, it indicates that the energy flux 
mostly comes from the inner region of the disc. 

The ratio of accretion power dissipated into the corona (/) 
as functions of the parameter Ae and the dimensionless BH spin 
parameter a* are shown in Fig. 2 and Fig. 3(a), and the magnetic 
pressure Pmag = ctoPtot is adopted in our calculations. From the 
above Fig.2 and Fig.3(a), we find that the global value of (/) is 
insensitive on the parameter Ae for Ae > 1, but the global value 
increases with the increasing the spin parameter a, . The curves of 
(/) versus accretion rate rh for the deferent BH spin parameter a» 
are plotted in Fig3.(b). The dimensionless accretion rate rh is de- 
fined as rh = M/MEdd in the calculations. As shown in Fig3b 
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Figure 1. The value of Q/Qo varying with r / r ms for the different values 
of Ae : Ae = (solid lines), Ae = 0.1 (dashed lines), Ae = 1 (dotted 
lines); The dimensionless BH spin parameter a„ = 0.95 is assumed. 
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Figure 2. The global value of (/) as a function of the parameter Ae for 
three different values of a» : a* =0 (solid lines), a* =0.5 (dashed lines), 
a„ = 0.9 (dotted lines), rh = 0.1, a = 0.3 is adopted in the calculations. 



the global value (/) is independent of accretion rate rh for a given 
BH spin parameter a* if the magnetic pressure P mag = cepPt ot is 
adopte d. We find that this result is consistent with those in ICaol 
(2008), though the central black hole is Kerr black hole in our 
model, while a Schwarzschild black hole in lCaol (2008). 

We plot the curves of the global value (/) versus different 
parameters with P mag = a P ga s and P mag = a \/ 'PgasPtot in 
Fig 4 and Fig 5, respectively. 

Inspecting the Fig.4 and Fig.5 we find that our disc-corona 
model calculations with P ma g = ctoP ga s show the integration 
fraction of accretion power dissipated into the corona (/) ~ 
0.270 - 0.298 with the BH spin pa rameter a , ~ 0.01 - 0.998, 
and the model with P mag = ao-y/ PgasPtot shows that (/) ~ 
0.285 - 0.315 with the BH spin parameter a, ~ 0.01 - 0.998 at 
rh < 0.05, for different values of Ae. 

From Fig.4(a) we find that the values of (/) depend on the 
parameters a* and Ae when Ae < 1. (/) decreases (increases) 
as Ae changes from to 1, for a, < 0.4 ( a* > 0.4 ). In fact, 
the magnetic fields connecting the plunging region to the disc can 
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Figure 3. (a) The global value of (/) varies with a„ for different values of 
Ae: Ae = (solid lines), Ae = 0.1 (dashed lines), Ae = 1 (dotted lines). 
m = 0.1 is adopted in the calculations, (b) The global value of (/} varies 
with m for different values of a*: a* = 0(solid lines), a, = 0.5 (dashed 
lines), a* = 0.9 (dotted lines), Ae = 0.1 and a = 0.3 are adopted in the 
calculations. 



(b) 

Figure 4. (a) The global value of {/) varies with a* for different values of 
Ae: Ae = (solid lines), Ae = 0.1 (dashed lines), Ae = 1 (dotted lines). 
m = 0.1 is adopted in the calculations, (b) The global value of (/} varies 
with m for different values of a*: a* = 0(solid lines), a* = 0.5 (dashed 
lines), a* = 0.9 (dotted lines), Ae = 0.1 and a = 0.3 are adopted in the 
calculations. 



exert stresses on the inner edge of the disc and transfer energy from 
the plunging region to the disc-coronae system. So the values of 
Qcor and Q all increase with the increasing Ae. Q cor and Q are 
also positively correlated to BH spins. For a* < 0.4, the value 
of Qcor increases more slowly than that of Q with the increasing 
Ae. Thus the ratio of accretion power dissipated into the corona, 
(/), decreases with the increasing Ae. For a* = 0.4, we find that 
the change of the magnetic torque has no influence on the ratio of 
the power dissipated in the corona to the total power. This critical 
BH spin parameter is about a* = 0.78 in Fig.3(a), and is about 
a, — 0.6 in Fig. 5(a). 

In addition, it should be noted that the value of (/) is also 
insensitive on the Ae for Ae > 1, as P m ag = cioPgas or Pmag = 
ceo \J PgasPtot is adopted in our calculations. 

As shown in Fig.4(b) and Fig.5(b), the value (/} decreases 
monotonically with the increasing accretion rate for m ^ 0.05 
. (/) can reach ~ 0.3 at the lower accretion rates. So it seems 
that th e low luminosity BH can b e fitted by our disc-corona model 
nicely. iMerloni & Fabianl d2002h also proposed that the magnetic 



corona should be stronger at low accretion rates, and their strength 
depends upon the nature of magnetic dissipation inside the disc. 



4 SIMULATION OF DISC SPECTRUM 

In the disc-corona scenario, the comptonized spectrum has been 
computed by some authors using different approaches in previ- 
ous works. Two kinds of approaches are used in order to get 
the emergent spectrum. One of common approaches is to solve 
the radiative transfer equat i on either numerically or an alytically 
dSunvaev & Titarchuklll980l : |Poutanen & Svenssonll996l). Another 
kind of approach is the Monte Carlo simula tion! Pozdny akovetajj 
1 19831; iGorecki & Wilczewskill 19841: Istern et al. Ill995l : lHualll997t 
iMa et alj|2006h . Recently. iMa et alj fcOOdT got the output spec- 
tra in the cases with and without the magnetic coupling effects ( 
Li 2000, Li & Paczynshi 2000, Li 2002a, Li 2002b ) by using the 
Monte Carlo simulation. 

In this paper, the steps of our simulations are: (i) sample a 
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Figure 5. (a) The global value of {/)varies with a» for different values of 
Ae: Ae = (solid lines), Ae = 0.1 (dashed lines), Ae = 1 (dotted lines). 
fn = 0.1 is adopted in the calculations, (b) The global value of (/} varies 
with m for different values of a* : a* = (solid lines), a* = 0.5 (dashed 
lines), a* = 0.9 (dotted lines), Ae = 0.1 and a = 0.3 are adopted in the 
calculations. 



seed photon from the cold disc; (ii) draw a value for its free path 
and test whether it can leave the disc or corona; (iii) simulate the 
interaction of the photon with the medium; (iv) repeat steps (ii), 
(iii) till the photon leaves the system of the corona and disc. 

In our disc-corona scenario, the dissipated power in the unit 
area of the disc, F is related to Q by F = Q — Q C or, and according 
to Stefan-Boltzmann law we have the local effective temperature on 
the disc expressed by 



T d (r) = (F/a SB ) 1/4 , 



(19) 



where <jsb is the Stefan-Boltzmann constant. The local radiation 
spectrum is defined by the Planck function: 



fl„(r) 



21) 



c 2 exp[hu/kT d (r)] - 1' 
thus the multicolor black-body spectrum of the disc is 



B v (r)2irrdr. 



(20) 



(21) 



The probability density of the seed photons can be written as 

P(r,E) = P(r)P(E), (22) 

where P(r) = L(r)dr/L is the probability of a photon emitted in 
the ring r ~ r + dr, and L is the total luminosity of the disc. P(E) 
is the nu mber density of photons having an energy E = hv that is 
given bv lPozdnvakov et alj ( 1983) 



P(E) = 



2C(3) 



b 3 E 2 (e 
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(23) 



where b = 1/Td , C(3) = 1.202 is the Riemann Zeta function. 

In our simulations,if the seed photon is scattered by the elec- 
tron in the corona, the optical depth that the photon travels between 
the i-th and (i+l)-th scatterings can be drawn with n — — In A, 
where ^ A sj 1 is a random number corresponding to a random 
event. The free path of photon can be drawn with 



— In XcttHc 



(24) 



where n e , a and <tt = 6.65 x 10~ 24 cm 2 are the number density of 
the electrons, cross-sections of scattering and Thomson scattering, 
respectively. H c and r c are the vertical height and optical depth of 
the corona. Since the electrons in the hot corona are relativistic, the 
cross-section of scattering a depends not only on the energy of the 
photon but also on the energy and direction of the electron. So we 
us e the cross- section averaged over the distribution of the electrons 
in lHual dl997l) to draw the free path of the photon. 

In the disc-corona system, hard X-ray photons irradiating the 
disc from the corona can be absorbed by the atoms in the disc, 
as well as being scattered by the free electrons. In this case, we 
choose the lesser of the two free paths that are drawn with a a and 
cr s to draw the free path of the photon, where a a and a s are re- 
spectively the cross sections of absoiptio n and scattering. The cross 
section s of absorption a a is taken from iMorrison & McCammonl 
dl983l) , and the cross sections of scattering can be given by the 
Klein-Nishina formula 



3<tt 
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,(25) 



where e is the energy of the incident photon in unit of electron-rest- 
energy. 

As the free path is known, the position that the photon arrives 
before the next interaction can be calculated. If the photon is out- 
side the disc-corona system, it will escape away from the system 
and its energy and direction are recorded. If the photon transfers 
from the disc to the corona or inverse, the point where the trajec- 
tory of the photon cross the interface between the disc and corona 
will be regarded as the next initial position of the photon in calcu- 
lations. 

In the simulating the interaction of photon and electron, we 
can sample a electron from thermal distribution, then calculate the 
energy and direction of scattered photon. In the calculation, we fol- 
low the calculation procedure described in P83.The bound-free ab- 
sorption of hard X-rays by the atoms in the disk will lead to ioniza- 
tion and vacancy, and induce emission of fluorescence lines with 
the probability called fluorescence yield, Y . In our simulations, if 
the photon is absorbed by a certain atom or ion, we can draw a 
random number A ( ^ A ^ 1) and compare it with the corre- 
sponding fluorescence yield. If A ^ Y, an emission line is brought 
out, whose direction can be sampled from the isotropic distribution. 
If A > Y, the photon vanishes and its trajectory ends. We only con- 
sider the 6.4 keV Fe K a fluorescence line in our simulations, due 
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Figure 6. The emerged spectrum from the disc-corona system. The total 
emissive spectrum and its black body,power-law and reflected components 
are shown in solid , dot-dashed, dotted and dashed lines, respectively. The 
radii of the inner and outer edges of the corona are taken as r^ n = r ms and 
Tout = 100r ms , and a* = 0.998, Ae = 0.1 and rh = 0.05 are adopted 
in the calculations. 



to the large abundance and cross-section of absorption of iron. If 
the energy of the absorbed photon is less than the iron K-shell ab- 
sorption edge, i.e. E < 7.1 keV, no emission line is produc ed, while 
Fe K a lines emanate with Y = 0.34 for E > 7.1 keV dMa et al.l 
120061) . 

In our stimulations, the half height of the disc is assumed to 
be H = 0.05r ms . All the spectra in this paper are given in forms 
of curves rather than histogram to make the curves look smooth. 
Fortran 95 is used in the calculations and Mathematica 7.0 software 
is adopted to plot the figures. 

The typical results for the spectra of our Monte Carlo simula- 
tions are given in Fig. 6. The spectrum of the disc-corona system 
comprises three components, the black-body spectrum formed by 
unscattered photons, the power-law spectrum formed by photons 
that escape from the corona after several times of inverse Compton 
scatterings, and the reflected spectrum characterized by the iron flu- 
orescence line and reflection hump. 

Spectra with different coronal geometry parameters are given 
in Fig. 7. We stimulate spectra with different heights of corona: 
He = 0.5r ms , He = r ms , He = 2r ms , in Fig. 7(a). It is found 
that the spectra steepen in the 2-200 keV range with the increas- 
ing height of corona. This result can be explained as follows: when 
the corona becomes thicker, the amount of photon-electron scat- 
tering increases, and the energy of the escaped photons increases. 
Therefore, hardening of the spectrum is anticipated. We present 
the spectra with different radii of the outer edge of the corona : 
Tout = 100r ms , r out = 50r ms , r out = 20r m3 in Fig.7(b). The 
fluxes increase with r ou t, since a larger disc can lead to more en- 
ergy from the coronal surface. The spectral profile changes slightly 
in this case. 

Spectra with different accretion rates rh, Ae and BH spin pa- 
rameters a„ are given in Fig. 8. From Fig. 8(a) we find that the 
fluxes increase as rh changes from 0.01 to 0.1, and the spectral 
profile changes obviously. As the accretion rate rh increases, the 
total gravitational energy dissipated is enhanced remarkably. In 
this case, more seed photons escape from the surface of the disc; 
these photons also have higher energies. It is found in Fig.8(b) 
that the fluxes increase as Ae increases from 0.1 to 2, since the 
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Figure 7. The simulated spectra of the disc-corona system, (a) Solid, 
dashed, and dotted lines correspond to the spectra with different heights 
of corona : He = 2r ms , He = r ms , He = 0.5r ms , for the radii 
of the inner and outer edges of the corona are taken as T^ n — Tms 
and r ou t = 100r ms . (b) Solid, dashed, and dotted lines correspond 
to the spectra with different the radii of the outer edges of the corona : 
r ut = 100r ms , r out = 50r ms , r out = 20r ms , for r in = r ma and 
He = r mB . as, = 0.5, Ae = 0.1 and fn = 0.05 are adopted in the 
calculations. 



magnetic torque can transfer energy from the plunging region to 
the disc-corona system. The energy distribution of the seed pho- 
tons, hence the emerged spectrum, also changes in this case. From 
Fig.8(c) we see that the fluxes increase with a„. Indeed, as a„ in- 
creases the inner edge of disc approaches the BH horizon rg = 
M(l + yl — ai) , and more gravitational energy can be released. 

M = 1OM , and P mag = ot P ga3 P to t is adopted through the 
calculations. 



5 DISSCUSSION 

In this work, we propose a disc-corona model, in which magnetic 
fields exert a torque on the inner edge of the accretion disc, and 
part of gravitational energy is dissipated in the hot corona. The to- 
tal gravitational power Q is derived from the thin-disc dynamics 
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Figure 8. The simulated spectra of the disc-corona system, (a) Solid, 
dashed, and dotted lines correspond to the spectra with different accretion 
rates : m = 0.1, m = 0.05, m = 0.01, for a, = 0.5, Ae = 0.1; (b) 
Solid, dashed, and dotted lines correspond to the spectra with different Ae : 
Ae = 0.1, Ae = 1, Ae = 2, for a, = 0.5, m = 0.05 ; (c) Solid, dashed, 
and dotted lines correspond to the spectra with different spin parameters : 
a, = 0.998, a, = 0.5, a* = 0.1, form = 0.05, Ae = 0.1. H c = r ms , 
fin = Tms and Tout = 100r ms are adopted in the calculations. 



equations and the global solutions are obtained numerically. It is 
found that the fraction of the power dissipated into the corona in 
the total for such disc-corona system increases with the increasing 
the dimensionless black hole spin parameter a«, but is insensitive 
on the parameter Ae for Ae > 1. 

We simulate the emerged spectra from the disc-corona sys- 
tem for the different parameters using the Monte-Carlo method. 
It is found that the spectral profile changes obviously with vary- 
ing the height of the corona and the accretion rate of the disc. The 
reasonable geometry of the corona is important for simulating the 
emerged spectra using the Monte-Carlo method. We adopt a slab 
corona in this paper. However the geom etry of the corona is still 
matter of debate, e.g. lStern et al. ]|l99i proposed a patchy corona 
that made of a number of separate active regions. Furthermore, the 
observe d correlation be t ween the photon index and the reflection 
strength Zdziars ki et al.l (1 19991) has demonstrated the need for fur- 
ther geometrical/dynamical parameters, such as the relativis tic bulk 
motion velocity of the coronal material (Ma lzac et al.ll200lb . 

Our model needs to be improved in other aspects. For exam- 
ple, the gravitational effects on the trajectories of photons need to 
be taken into account, the cooling of synchrotron radiation should 
be considered and the ray-tracing should be used in our calcula- 
tions. 

The magnetic pressure P mag plays an important role in our 
model. In fact we can get the magnetic energy density from the 
coronal power Q cor as the coronal geometry is assumed. Then the 
small-scale magnetic field Qdynamo in the corona can be derived 
from the magnetic energy density. 

On the other hand, models and simulations of jet produc- 
tion ilBlandford & Znaieklll977h iBlandford & PavneHl982t iMeied 
ll999l) show that it is the poloidal component of the large-scale mag- 
netic field which mainly drives the production of powerful jets. 
Several theoretical models have been proposed for acceleration and 
collimation of jets, which can be divided into two main regimes, the 
Poynting flux regime and the hydromagnetic regime. Both regimes 
are related to a poloidal magnetic field threading the disc or BH, 
from which energy and angular momentum are extracted. In the 
Poynting flux regime, energy is extracted in Poynting flux (i.e. 
purely electromagnetic energy), but in the form of magnetically 
driven material winds in the latter regime. Furthermore, some au- 
thors have agreed that jet formation should inv olve an accretion 
disc threade d by a large-scale magnetic field dLivio et all 1 19991 : 
lMeiej|l999h . 

Though the origin of large-scale magnetic field is still 
under controversy, so me previous works dTout & Pringle Il996l : 
iRomanova et al.f l998) have proposed that the large-scale field can 
be produced from the small-scale field created by dynamo pro- 
cesses. The length scale of the fields created by dynamo processes 
is of the order of the disc thickness H, and the poloidal component 
of the magnetic field is given approximately by 



Bp ~ (H/r)B. 



dynamo j 



(26) 



If the field is created in the thin accretion discs (H <C r), the large- 
scale field is very weak. For the ADAF cases, the disc thickness 
H ~ r and the poloidal component of the magnetic field shall be 
stronger. In our disc-corona scenario, the energetically dominant 
corona are the ideal sites for launching the powerful jets/outflows 
dMerloni & Fabianll2002l) . The large-scale magnetic fields created 
by dynamo processes in the corona are significantly stronger than 
the thin disc due to the corona being much thicker than the cold, 
thin disc. So the corona can power a stronger jet than the thin disc. 
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We shall discuss the acceleration of jet in the Poynting flux regime 
and the hydromagnetic regime in future work. 
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